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Abstract 

We define a family of Schrodinger Functional renormalization schemes for the four-quark 
multiplicatively renormalizable operators of the AF = 1 and AF = 2 effective weak 
Hamiltonians. Using the lattice regularization with quenched Wilson quarks, we compute 
non-perturbatively the renormalization group running of these operators in the continuum 
limit in a large range of renormalization scales. Continuum limit extrapolations are well 
controlled thanks to the implementation of two fermionic actions (Wilson and Clover) . The 
ratio of the renormalization group invariant operator to its renormalized counterpart at a 
low energy scale, as well as the renormalization constant at this scale, is obtained for all 
schemes. 



1 Introduction 



In the quest for accurate quantitative predictions in the non-perturbative sector 
of the Standard Model, non-perturbative renormalization has become an essential 
element of lattice QCD calculations. Two schemes are currently in use in many 
applications, namely the RI/MOM and the Schrodinger Functional (SF) |2j. In 
the latter scheme the scale evolution of (matrix elements of) renormalized opera- 
tors can be traced non-perturbatively over a wide range of scales. The validity of 
perturbation theory at high scales can thus be verified and one may then convert 
perturbatively to one of the commonly used continuum schemes, such as the MS 
scheme of dimensional regularization. Alternatively, one may use low order per- 
turbation theory to extrapolate from high to infinite energies, where the so-called 
renormalization-group-invariant (RGI) operators are defined. In any case, pertur- 
bation theory is only used in the high energy regime where it may be safely applied. 

These techniques have been applied to study the scale evolution of various 
physical quantities, such as the QCD gauge coupling, the quark mass HJ 00 E] 
and the moments of pion or nucleon structure functions 7 j (both in the quenched 
approximation and for two dynamical flavours), as well as matrix elements of the 
heavy-light axial current, with heavy quarks treated in the static approximation [S]. 
The present work is a first step towards the extension of this SF renormalization and 
renormalization group (RG) evolution programme to four-quark operators relevant 
for weak matrix elements. These arise in the OPE as the low energy QCD contri- 
bution in weak interaction transitions. They are key elements for the determination 
of the CKM unitarity triangle (and the subsequent understanding of CP-violation). 

We specifically investigate the renormalization of two dimension-six operators 
with a "left-left" Dirac structure and four fermions with distinct flavours: 

±{MxhflMx)) Wzfr^MaO)] (L1) 
= o ± -o ± 

W VV+AA ^VA+AV ' 

where 7^ = 7^(1 — 75). The last expression implicitly defines the parity-even and 
-odd components of O hL , in fairly standard notation. In a chirally symmetric reg- 
ularization, the operators 0^ L are multiplicatively renormalizable. As we opt for 
lattice regularizations with Wilson fermions, the loss of chiral symmetry generates 
extra parity-even counterterms with finite mixing coefficients. The parity-odd com- 
ponents Ova+av are protected against the generation of parity-odd counterterms 
by discrete symmetries. For a full account of these renormalization properties see 
refs. [H1E3J- I n the present work we focus on these multiplicatively renormalizable, 
parity-odd operators Ova+av 

Once the four generic flavours are identified with specific physical flavours, 
the corresponding weak matrix elements give rise to a variety of phenomenology. 
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For example, identifying ip\ and ^3 with the strange (or bottom) quark and ip2 
and ?/>4 with the down quark, we obtain the operator O^a+av mediating AF = 2 
transitions (K° — K° and B° — B° oscillations) in the tmQCD lattice regularization 
framework jllj . If ipi is a strange quark field and the others are suitably chosen 
light and charmed quarks, we are looking into AS = 1 operators O^a+av mediating 
K —* 7T7T transitions. Our results completely determine the renormalization of the 
operators mediating the AI = 3/2 channel, whereas for the AI = 1/2 transitions 
only logarithmic divergences are removed. This multiplicative renormalization is 
sufficient in the limit of SU(4) flavour symmetry, of which the chiral limit is a special 
case. Upon explicit breaking of this symmetry by the masses of the heavier quark 
flavours, the renormalization programme of the AI =1/2 channel is only complete 
once the question of mixing with operators of lower dimension has been addressed. 
This mixing is beyond the scope of the current work. 

The paper is organised as follows: In sect. 12 we present a general discussion on 
the RG running of correlation functions of composite operators, leading to the defini- 
tion of the corresponding renormalization group invariant (RGI) operators. In sect. 01 
we introduce the SF renormalization schemes used for the four-fermion operators in 
question, define the operator step scaling functions (SSF), discuss their properties 
and show how the operator RG running can be obtained non-perturbatively from 
them. In sect. |l]we present our non-perturbative computation of the SSF. Our re- 
sults have been obtained in the quenched approximation. We used both the standard 
Wilson quark action and its 0(a) improved version, with the Sheikholeslami-Wohlert 
or clover term ^2] (henceforth referred to as Clover action). Once extrapolated to the 
continuum limit, the SSF is used in order to obtain the ratio of the RGI operator to 
its renormalized counterpart at a hadronic low energy scale. In sect. 03 we compute 
the operator renormalization constants at this hadronic matching scale. Finally in 
sect. we discuss our conclusions. Some technical points have been relegated to 
appendices. Preliminary results had already appeared in refs. |13j . 

In a companion paper ^1] the same calculation has been performed in perturba- 
tion theory. The lattice SF schemes have been matched to a standard MS continuum 
scheme, at 1-loop in perturbation theory. Combined with the known NLO results of 
the operator anomalous dimension in the continuum reference scheme, these results 
give the NLO estimate of the operator step scaling function. These in turn are used 
in the present work, since part of the calculation involves the RG running of the 
operators at very high scales in the SF schemes, where NLO perturbation theory 
may be safely applied. 

2 Callan-Symanzik equations and RGI operators 

Our starting point is the Callan-Symanzik equation expressing the RG running of 
correlation functions under a change of renormalization scale \x. Our exposition and 
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notation follows closely that of Refs. [3J [EJ. We first consider an arbitrary bare 
n-point correlation function 



G(x 1 , . . . ,x n ;g ,m ,i) = (Oi(si) ■ ■■O n (x n )) 



r-l 



V[i>^]V[U]e- s 1 {x 1 )---O n {x n ) 



2.1) 



where the Oi are local gauge invariant composite operators and Z is the QCD 
partition function. A regularization such as the lattice (with ultraviolet cutoff a -1 ) 
is implied. For simplicity we assume that all space-time points are separated; i.e. 
Xi 7^ xj for i ^ j; i, j = 1, . . . , n. The dependence on the bare parameters go,mo^ of 
the theory has been indicated explicitly. As this is quite cumbersome, occasionally 
we will omit some of the arguments, in order to simplify the notation. The subscript 
f = 1, . . . , Nf borne by the mass indicates flavour (the mass matrix will be henceforth 
assumed to be diagonal). 

It is adequate for the purposes of the present work to consider only multiplica- 
tively renormalized operators (i.e. no operator mixing occurs). Their renormalized 
correlation functions at scale u can be written as 



G R (x 



Y[z Ol {g ,an) 



.1=1 



G(xi,. . . ,x n ;g ,m ,{) ■ (2.2) 



We denote the renormalized coupling by or and the renormalized quark masses by 
m^f. The operator renormalization constants Zo t are determined by imposing n 
renormalization conditions on suitably chosen correlation functions of the operators 
Oi(x), at scale /i. In the context of the present discussion these conditions need not 
be specified; it is crucial however to keep in mind that they are imposed in the chiral 
limit i.e. we are only considering mass independent renormalization schemes. 

The correlation function defined in eq. (|2.2|) fulfills a Callan-Symanzik equation 
which determines the RG running of the operator in question: 



5 Nf d 

+ P(ffR)-£Z7 + r (f r) mR < f 
f=l 



dfj, 



dg 



R 



dm 



K.f 



G R = O(a) 



(2.3) 



where /3(or) in the Callan-Symanzik function, t(or) the quark mass anomalous 
dimension and 7o j (flR) the anomalous dimension of operator Oi, which is related to 
its renormalization constant through 



7o(ff(M)) 



lim 



fi^-Zo{9o,afjL) 



Zo(go,a(-i)~ 



(2.4) 



Since the renormalization scheme we are working in is mass independent, /?, r and 
7o depend only on the renormalized coupling. Their asymptotic expansions at small 
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values of the coupling are given by 

0(g) 9 ~° ~g 3 (bo + hg 2 + b 2 g A + ...) , (2.5) 

r{g) -g 2 (do + d l9 2 + d 2 g A + . . .) , (2.6) 

I \ 9^° 2 ( (0) . (1) 2 , (2) 4 , \ /r, v \ 

70(5) 5 (7o +To 9 +r g +■■■) ■ ( 2 - 7 ) 



Running parameters are then defined as usual: 

q % = mq)) > q l^~ = T{v{q)) mdq) ' (2 ' 8) 

supplemented by the boundary conditions 

<7(a0 = 3R> mf(^) = m R,f • (2.9) 

To define RGI composite operators, we start with the formal integration of 
eq. (|2.3|) . yielding 

G R (xi, . . .,x n ;n',g(ji'), m(n')) = 

n (2 10) 

where Ui is the evolution function 



This function describes the RG evolution in the continuum limit of the renormal- 
ized operator (Oj)r between the renormalization point and an arbitrary scale //, 
namely: 

(Oi) R (x; m') = /i) (Oi) R (x; //) . (2.12) 
It can easily be seen to satisfy the RG equation 

q 9W^ =io{mmq ^ ) (213) 

with initial condition 

Ui(fi,fj,) = l . (2.14) 

From eqs. (|2.2[1 and eq. ()2.10|) we can also express it as a ratio of renormalization 
constants 

tt 1 1 \ v z Oi(go,afi r ) 
°- >0 z o Ago, aii) 
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The RGI operator could in principle be obtained by splitting the r.h.s. of 
eq. (|2.11j) into two integrals (one from g(n') to g = and one from g = to g(n)) 
and subsequently bringing the //-dependent integral on the l.h.s. of eq. 
problem is that Ui(fi' , /x) diverges logarithmically in the limit \J 



oo. 



The 
This is 

most clearly seen upon considering the asymptotic expansions admitted by (5 (cf. 
eq. (|2.5jl ) and 70 (cf. eq. (j2.7|0 at small values of the coupling. Hence, we proceed 
by casting eq. (|2.12|) in the form 



IV) 

47T 



-7^/(2feo) 



x exp 



em 



-7^/(2feo) 



dg 



gin 1 ) 



7o(g) 



(o). 
To 



P(g) b g 



Or(x;h) 



(2.16) 



What has been achieved is the finiteness of the r.h.s. of eq. (|2.16|) as // — ► 00 (i.e. 
g(n') — > 0) for any value of fj,. Moreover, since there is no ^-dependence on the l.h.s., 
also the r.h.s. is /U-independent. Hence, taking the limit // — > 00 of eq. (|2.16|) . we 
define a RGI quantity as 

d(x) = Z (fJi)0 R (x;tJ,) , (2.17) 

where we have introduced 



Zo(fJ-) 



47T 



It must be stressed that the RGI operator 0(x) defined above is (unlike Or(x; /i)) 
independent of the renormalization scheme and scale. The expressions (|2.17|) and 
(|2.18|) for O(x) are an exact result, in close analogy to the ones reported in ref. [1] 
for the RGI quark mass 

M, = W M-V)r^W{-£ £)} , ( , 19) 

and the RGI scale 

A = /i(6o5 2 (/^))" bl/(2bi) exp 

f rf&O ^ ^ n ^ (2 - 2U) 

x exp ^ - / dfif 








Notice that the only arbitrariness in the definition of these RGI quantities is a 
constant overall normalization factor. For composite operators in Eq. ()2.17|) we 
have adopted the normalization usually employed in the definition of the RGI kaon 
mixing parameter Bk- 
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3 RG running of four-fermion operators 



Having exposed the general principles for the RG behaviour of multiplicatively renor- 
malizable composite operators, we now pass to the specific case of interest. This 
refers to the dimension-six composite operators of four distinct quark flavours 

Ova+avO) =\{ [(i>nM (^37^75^4) + (^17^75^2) (^37^4)] ^ 3 1 , 

± [i^ll^A) (^37^75^2) + (^17^75^4) {^1^2)] } , 

which are known to be multiplicatively renormalizable E3] • In this section we 
define the correlation functions of interest, the renormalization conditions imposed 
and the step scaling functions of the operators Ova+av 

As anticipated, we opt for the lattice Schrodinger functional (SF) formalism |17l 
I18[ I19|. We regularize QCD on a lattice of extension L 3 xT (here T = L always) with 
periodic boundary conditions in the space directions (up to a phase 9 for the fermion 
fields) and Dirichlet boundary conditions in the Euclidean time direction |18[ 119] . 
Otherwise the lattice gauge and fermionic field actions are of the standard Wilson 
type; the Clover 0(a) improved version of the fermionic action is also used. The 
operators Ova+av i x ) are defined locally on the lattice; i.e. all quark fields live at the 
point x. 

The paper follows closely the notation of ref. (201) to which the reader is referred 
for unexplained notation. 

3.1 Schrodinger Functional correlation functions 

Bare composite operators are defined at both time boundaries in terms of the bound- 
ary fields C and (' of refs. QH HOI EH, 

o 12 [r] = a 6 ^c"i(x)rc 2 (y), 

x,y 

o'i 2 [r] = a 6 ^C((x)rc 2 (y). 3 ' 2 

The indices 1,2 label distinct flavours. Unprimed fields are defined on the xq = 
boundary, primed ones on the xq = T one. There are two allowed independent 
choices for the Dirac matrices, namely T = 75 and T = 7^ (with k = 1,2,3). This 
is due to the SF Dirichlet boundary conditions of the quark fields, which involve 
positive and negative projection operators P± = |(1 ± 70) HJ. The presence of 
these projectors implies that boundary sources with other Dirac matrices T either 
vanish or are identical to those with 75 or 7^. 
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The bare correlation functions of the four-fermion operators are now chosen as 
follows: 



^A.iB.rc]^) = ^3-(02i[r A ]045[r B ] Ot A+AV (x) 0' 53 [T C ]) . (3.3) 

These are depicted, in terms of valence quark lines, in Fig. ^ Since the boundary 
operators defined in eqs. (|3.2|) involve sums over all 3-space at each time boundary, 
translational invariance implies that the above correlation functions depend only 
on time. A few words are in place in order to motivate the choice of this rather 
complicated quantity, involving three composite operators at the boundary. As 
stated above, the boundary operators can either be Ci2[75] or Oi2[7fc] (and C 12 [75], 
0' l2 [lk\)- Moreover, the bulk operators Ova+av( x ) are parity-odd. These facts, 
combined with the requirement that cubic symmetry be respected by the correlation 
functions, give as simplest possibilities the following five, in principle independent 
correlation functions: 



F 1 ± (xo)=^ 5i75i75] (x ), 



7± 

1 3 

j,k,l=l 
1 3 



3 [75,7fc,7fcP v/ 7 (3_4) 

1 3 

k=l 
1 3 

5 V yJJ o [7fc,7fc,75] v u; 



k=l 

We will also need the boundary-to-boundary correlation functions 

fx = -^(OM om), 

3 (3.5) 

k=l 

In terms of valence quark propagators, these correlation functions are depicted in 
Fig. El 

In practical simulations, these correlation functions are computed as traces 
of the boundary-to-bulk valence quark propagators Hf(x) and H'^x), defined in 
ref. [22]. 
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Figure 1: Four-fermion correlation functions in a finite spacetime volume with 
Schrodinger Functional boundary conditions. The lines denote valence quark prop- 
agation of different flavours; the double dots in the bulk denote the four-fermion 
operator. 



3.2 Schrodinger Functional renormalization schemes 

Before turning to the renormalization of four-fermion operators, we recall that SF 
renormalization schemes are mass independent; i.e. renormalization is performed 
in the chiral limit. The renormalization scale is set at fi = 1/L; the renormalized 
coupling g(l/L) (defined in ref. [Hj) and quark mass m(l/L) (defined in refs. |4*|IT5]) 
are then only functions of this scale. 

Upon removing the ultraviolet cutoff (i.e. a — » 0), the bare correlation functions 
defined in eqs. (|3.3U3.5|) diverge logarithmically. The divergence due to the boundary 
fields is removed by considering suitable ratios of correlation functions. Several 
choices can be made, giving rise to different correlator ratios. In the present work 
we will be considering the following nine specific cases 

\ («o) = — §75— 1 = 1, ... ,5 , 
fi 

h efa) = %^ , (3.6) 

, ± / \ F±(x ) . 
^1+4(^0) = —Yf 2 — 1 = 3, 4, 5 , 
fi fa 

which renormalize as the four-fermion operators Ova+av themselves: 

h% s (x ;ij) = Z^ A+AV . s (go,a/J,)ht(x ;g ) s = l,...,9. (3.7) 

The above renormalization constants are fixed by imposing the following renormal- 
ization conditions on the correlator hf^ on time-slice xq = L/2 (for all s = 1, . . . , 9) 
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Figure 2: Boundary-to-boundary correlation function in a finite spacetime volume 
with Schrodinger Functional boundary conditions. The lines denote valence quark 
propagation of different flavours. 

at scale \x = \jL and fixed renormalized coupling g 2 (l/L) = u in the chiral limit: 



■^VA+AV;s (g ,an)hf(x ;go) 



hf{x ;g ) 



(3.8) 



90=0 



i.e. at tree level Zy A+Av . s = 1 by construction. We will always impose eq. 1)3.8)1 at 
9 = 0.5 [HE]- The nine correlator ratios chosen above give rise to nine in principle 
distinct SF renormalization schemes for each of the two operators. 1 

The above construction hinges on a theory of five flavours, the fifth one being 
a spectator quark. This, however, is only apparent. As already mentioned in the 
introduction, we make contact with a specific weak matrix element by judiciously 
attributing specific physical flavour labels to the above five nominal flavours. For 
example, the identifications 



4>l = 03 



02 = 04 



d. 



^5 



II, 



(3.9) 



lead (up to an irrelevant factor of 2 arising from a doubling of Wick contractions) 
to the renormalization of the "left-left" operator Ova+avi which mediates A 5 = 2 
transitions. 



3.3 Step scaling functions 

The step scaling functions (SSF) of the four-fermion operators of eq. 1)3. lj) are defined 

as 



VA+AV; 



,(g ,a/2L) 



■^VA+AViB 

(g ,a/L) 



(3.10) 



m=0, g 2 (l/L)=u 



1 Some considerations concerning the independence of the different schemes can be found in Ap- 
pendix^ 
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This is in close analogy to the quark mass case [3]; i.e. Eva+av ;s is defined in the chiral 
limit m(go) = 0, for a lattice of a given resolution L/a and at fixed renormalized 
coupling g 2 (l/L) = u. The precise definition of the current quark mass m(go) can 
be found in ref. [I]. The lattice SSF Sv A+AV;s is not unique: it depends on the details 
of the lattice regularization (e.g. the type of lattice action chosen, the level of 0(a) 
improvement etc.). It has, however, a well defined continuum limit, which should 
be unique (i.e. universality should hold). We denote the continuum SSF by 

^VA+AV;» = J™ S VA+AV ;S ("> a / L ) • ( 3 - U ) 

In terms of the operators' evolution function J7 V a+av- s and anomalous dimension 
7va+av; S the SSF can be written as (cf. eqs. IJ2.11I2.15J) ) 

± / \ rri ( 1 l \ \ f~ gil/2L) 7vA + AV: S (g) A 1 

<Wv» = U VA+AV;s (-, j) = exp | m d,| . (3.12) 

Thus the physical meaning of o"v A+AV;s emerges readily from the above as the operator 
evolution function between two scales differing by a factor of 2. It is a quantity closely 
related to the anomalous dimension of the corresponding operator. 

We stress that the operator anomalous dimension is scheme dependent. Its 
perturbative expansion is known to two-loop order; the universal one-loop coefficient 
is 

± 1 6(iVTl) (3 13) 

7vA+AV-±^ 2 N , (A-WJ 

(N being the number of colours) while the two-loop coefficients 7va+av ;s have been 
calculated in 14 for the schemes defined by eqs. ()3.8() . 

Finally, from eq. (|2.18|) we immediately obtain the following expression for the 
factor relating the renormalized operator at a scale [i with its RGI counterpart: 



^VA+AV;s(^) 



5 2 (/i)l" 7 vA+Av/( 2fe o) 



47T 

x exp <{ - / dg 



x 

±(o) ( 3 ' 14 ) 

' 7vA+AV;s (g) 7vA+AV;s 







It is the aim of the present work to provide accurate estimates of the above quantity 
in all nine schemes and for a large range of scales fx = 1/L (albeit in the quenched 
approximation) . 

It is useful to keep in mind that the only flavour dependence of ^va+av-s(/ x ) i s 
through Nj; there is no dependence on the values of the physical quark masses, as 
this quantity in defined (and computed) in the chiral limit. Thus it can be readily 
used in the renormalization of various physical matrix elements. 



10 



3.4 RG running of four-fermion operators 

Once the step scaling functions <7va+av ;s nave been computed through numerical 
simulation, the ratio of renormalized correlation functions involving the four-fermion 
operators of interest between the minimum and maximum renormalization scales 
covered by these simulations can be worked out. In order to be consistent with the 
notation of ref. [I], we denote the former by // m i n = (2L max )~ 1 . The ratio in question 
is then obtained in two steps: 

First the SSF of the gauge coupling 

a(u) = g 2 (l/2L)\. 2{1/L)=u , (3.15) 

computed in pHH], is used in order to determine the correspondence between renor- 
malized couplings and renormalization scales. This is done through the recursion 

ui = cr(u/ + i) , (3.16) 

with uq = <? 2 (l/£ m ax) = 3.48 the initial value. 2 

Second the SSF Cva+av ; s; known non-perturbatively, is used for this sequence of 
couplings in order to compute the quantity 

fe-i 

^VA+AV;s(Mmm) 2 // m in) = I I [^Va+AVis («,)] • (3.17) 

1=0 

The number of recursion steps h has to be chosen so that the (large) scale 2 /imin 

lies 

in the range covered by the computation of the SSF. In practice (cf. section 0J), it is 
safe to take k = 7, which means that 2 fc /x m ; n is deep in the region where perturbation 
theory can be expected to apply. 

The final step in our calculation is the computation of the RG running factor 
of eq. (J3.14|) at the renormalization point fj, m i n , written as a product: 

■^VA+AV;s(^ min ) = ^VA+AV;s(2 /^miii) ^min)-^vA+AV;s(^ ^min) • (3.18) 

The first factor on the r.h.s. is known from eq. (|3.17() . The second factor, which 
involves a (presumably) perturbative scale 2 k ^ m - m , is calculated from eq. 1)3. 14|) with 
the NNLO and NLO perturbative expressions of /3(g) and 7va+av ; s(<50> respectively. 
Clearly the underlying assumption is that the truncation of the perturbative series at 
NLO is safe at this scale. This is a scheme dependent statement. The perturbative 
results of ref. |H] indicate that, for schemes s = 1, 3, 7, the NLO coefficient of 7va+av ;s 
is of the same sign and much smaller than the LO one; s = 1 is the scheme with 
the smallest NLO coefficient. In all other schemes the relative sign is negative. For 

2 This initial value uo = 3.48 corresponds to L max /ro = 0.738(16); the initial calculation was 
performed in ref. |23) while the above result is obtained in the more recent ref. |24j . 
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7va+av;s the ratio is negative for all nine schemes, with s = 8 the smallest in absolute 
value. Conservatively, we indicate s = 1 and s = 8 as the most suitable schemes 
for operators Ova+av and Ova+av respectively. We stress that these choices are only 
dictated by the behaviour of the NLO perturbative results; the non-perturbative 
computation is equally reliable for all nine schemes considered. In any case, we have 
carried out our computations for all schemes. 

4 Non-perturbative computation of the step scaling function 

In this section we present the computation of Ey A+AV;s and its extrapolation to the 
continuum limit. We also obtain estimates of the corresponding RGI quantity (or, 
more precisely, of the expression of eq. (|3.18j) at a hadronic scale). The method of 
computation parallels closely that of refs. [3J[S] for the SSF of the quark mass S P . 

4.1 Wilson and Clover actions 

We have used both the standard Wilson action and its 0(a) improved version 
(Clover) in our simulations. Our notation is fairly standard; = 6/<7q is the inverse 
coupling and k = l/[2amo + 8] is the hopping parameter. At fixed bare coupling 
we define k ct as the value where the PCAC quark mass m(go) of ref. 0] vanishes. 
Following [3], the computation of done at 9 = 0. 

The 0(a) Symanzik improvement of the Schrodinger Functional has been worked 
out in refs. [2U1 [^J. For the pure gauge action, it amounts to modifying it 
by introducing time-boundary counterterms proportional to [ct(<7g) — 1]. For the 
fermionic action we must introduce the well-known clover counterterm in the lat- 
tice bulk, proportional to c sw (go), and time-boundary counterterms proportional to 
[ct(og) — 1]. Correlation functions of composite operators may then also be 0(a) 
improved by including in their lattice definition the appropriate higher dimension 
counterterms. Since for dimension-six operators, such as the ones of eq. ()3.1|) . there 
are several dimension-seven counterterms, we will not pursue operator improvement 
in this work. 3 

The improvement coefficient c sw has been computed non-perturbatively for a 
range of values of the bare coupling go ; see ref. |26| . The coefficients ct and 6t are 
known only in perturbation theory, to NLO |27j and LO 2X| respectively: 



In the present work we will distinguish two approaches to the continuum limit: 
3 0(a) improvement along the lines of ref. |25j is not readily applicable in the SF framework. 



ct( 5 o) = 1 - 0.089^ -0.030^ 
ct(g 2 ) = 1 - 0.018<7 2 . 



(4.1) 
(4.2) 
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(i) What we call "Wilson action results" (or "Wilson case" for short) consists in 
setting c sw = . Moreover, we set ct = 1, while the one- loop value 4 (eq. (|4,1|) 
truncated to 0(<?q)) is used for ct- 

(ii) What we call "Clover action results" (or "Clover case" for short) consists in 
using the Clover action with a non-perturbative c sw . The one- loop value from 
eqs. (|4.1[) and (|4.2|) is used for c t and c\ respectively. 

In both cases the four-fermion operator is left unimproved, so the dominant dis- 
cretisation effects are expected to be 0(a). We note, however, that the correlation 
functions (|3,4j) are 0(a) improved at tree- level, implying that all 0(a) counterterms 
to the local four-quark operators vanish at this order. Thus, for the Clover case we 
are left with discretisation errors which are 0{g$a). 

4.2 Continuum limit of the step scaling function 

For both the Wilson and Clover action the lattice SSF T^ A+AV . S have been evaluated 
at 14 values of the renormalized coupling g(l/L), each for four lattice resolutions 
L/a = 6,8,12 and 16. The tuning of (3 at the four L/a values, corresponding to 
a fixed renormalized coupling g 2 (l/L) = u, has been taken over from ref. [I]. The 
values of 

^cr taken from refs. HJ[71|!)|. The typical statistics accumulated for small 
lattices is of several hundred configurations. For the largest lattices the number of 
configurations ranges from around 60 at the weaker couplings to around 200 at the 
stronger ones. It has to be stressed that Wilson and Clover data have been obtained 
from independent ensembles of gauge configurations. 

A full collection of our raw data for Sv A+AV;s is available from the authors upon 
request. In Tables IC . 1I( L4I we present our results for Sy A+AV;1 and £y A+AV;8 ; see the 
discussion after eq. (|3,18|) for a motivation behind this choice. The quality of the data 
for the other schemes is comparable. The SSF £va+av ;s must be extrapolated to zero 
lattice spacing a/L (at fixed gauge coupling) in order to obtain its continuum limit 
counterpart fTvA +A v ;s - Since the four-fermion operators have not been improved, we 
expect the dominant discretisation effects to be O(a) both for the Wilson and Clover 
action data and thus a linear behaviour in a/L. Nevertheless we have performed fits 
on both datasets with two ansatze 

S va+av;sK a/L) = (Ty A+AV . s (u) + p{u){a/L) , (4.3) 
S va+av;s(^ a/L) = CTva+av;s(«) + p{u){a/L) 2 . (4.4) 

4 This is a choice of convenience: it is important to know for renormalization purposes (see 
eq. H3.10p below) the dependence of the Schrodinger functional renormalized coupling g(l/L) on 
the bare coupling go- This dependence is known non-perturbatively [H1IU for the pure Yang-Mills 
action with this c t value. In any case, the choice for c t has no bearing on the order of leading lattice 
artifacts. 
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An issue raised in refs. [UEj is the number of data points which should be included 
in each fit. In those works the L/a = 6 results were dropped from the fits, being too 
far from the continuum limit. We have performed fits with all data (4-point fits) 
and also without the L/a = 6 data (3-point fits). This means that we have applied 
a total of four fitting procedures (the two ansatze of eqs. Ij4.3|4.4|) . each for a 3- and 
a 4-point fit). 

The details related to the continuum limit extrapolation are presented in Ap- 
pendix El Prom that discussion we conclude that a conservative choice consists in 
performing 3-point fits (i.e. drop the data computed at the largest lattice spacing) 
which are linear in (a/L). The 1-loop perturbative discretisation errors have been 
divided out of £va+av ;s in the Wilson case. Moreover, following ref. [7j, we constrain 
the fits to the Clover and Wilson action data (at a given renormalized coupling) 
to have a unique continuum limit. 5 The outcome of this procedure is illustrated in 
Figs. IC.llC.2l for the two schemes of reference, and reported in Tables IC.5|C.6I for 
all schemes. We consider results obtained from these combined fits to be our best, 
and use them in the next step of the analysis. 6 



4.3 Continuum step scaling function and RG running 

The previous analysis has yielded accurate results for the continuum SSF Cva+av;s 
for a wide range of renormalized couplings. The data in this range of couplings can 
be represented by a polynomial of the form 

N 

± 



°va+av 



{u) = l + Y J stu n . (4.5) 



This ansatz is motivated by the form of the perturbative series. In perturbation 
theory the first two coefficients are known: 

±(°) 



sf =7va7av1h2, (4.6) 

4 = 7va ( ?av ; >2 + [^(TviSw) 2 + &o7vi?Av] ( ln2 ) 2 ' ( 47 ) 

The LO coefficient is universal, while the NLO one is scheme dependent and has 
been calculated in ref. |14j . As a result of the rather strong scheme dependence of the 
NLO anomalous dimension, also varies significantly between different schemes. 
In Figs. IC.3IC.4l (see left columns only) we compare the LO and NLO perturbative 
predictions for the SSF to the non-perturbative results of the present work. We 

5 This universality assumption has been thoroughly tested on our data for the SSF of the quark 
mass in 

6 We have also, in the spirit of ref. I28|. studied the impact of one-loop cutoff effects on the 
extrapolations. Some details are provided in Appendix 151 
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observe that, while the LO results are close to the non-perturbative ones at least for 
weak couplings, the NLO corrections show marked disagreement for certain schemes. 
This simply indicates poor convergence of the NLO perturbative series for some 
schemes. 

The values of the coefficients s^ of eq. 1)4.5)1 have been obtained through a 
suitable fitting procedure (see below for details). We can then compute the running 
of the composite operator between the scales // m i n and 2 fc /x m i n as explained in sect. 13.41 
(cf. eq. 1)3.17)) ). As input for the recursion in eq. 1)3.16)) we use the SSF of the 
renormalized coupling and its fit to a polynomial 



cr[u 



i + Yl a ^ 

n=l 



(4.8) 



as obtained in refs. H], with u\ and a% fixed from PT and 03, 04 kept as fit 
parameters. Then we apply eq. 1)3.17)1 with k = 7 iteration steps (corresponding to 
the range of scales covered by our simulation), and finally eq. ()3.18)) to obtain the 
RGI renormalization factor Z VA , A y. s (/x m i n ). The reliability of the computation of 
both factors on the r.h.s. of eq. ()3.18)) may in principle be compromised by a poor 
convergence of the perturbative series at NLO, which is indeed the case in some 
schemes. In particular the first factor could be affected if the coefficient S2 is kept 
fixed to its NLO value in the fit. The second factor could also clearly be affected, 
since it is calculated to NLO. 7 

Before addressing these issues, we discuss how we obtain a faithful fit to our 
data for Cva+av ;s ( u )) based on eq. 1)4.5)1 . We keep the first order coefficient si fixed 
to its perturbative value and perform a series of fits: 

(A) one-parameter fits with with S2 a free parameter; 

(B) two-parameter fits with S2 and S3 as free parameters; 

(C) one-parameter fits with S2 fixed from PT and S3 a free parameter; 

(D) three-parameter fits with S2, S3 and S4 as free parameters; 

(E) two-parameter fits with S2 fixed from PT and S3, S4 as free parameters. 

The results of these fits are summarised in Tables IC.7IC.8l We see that the SSF 
<7va+av ;s is well fit in all cases (x 2 /d.o.f. ~ 1). Also the SSF o"va+av ;s is always 
modelled well by the fitting curves, with the only exception of Fit C in schemes 1 
and 7, where the x 2 /d.o.i. is slightly higher. In any case, it appears that even in 
those schemes where the NLO RG running does not match the NP one, the fits are 
satisfactory, as the effect of the fixed NLO value of s^ is compensated by the higher 
order free parameters. 

7 In practice the calculation of this second factor is performed by numerically integrating the 
first of eqs. 12.81 (with the /3-function given at 3 loops) followed by numerical integration of the 
exponent in eq. I|3.14[l (with the operator anomalous dimension given at 2 loops and the /3-function 
at 3 loops). 
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The results for the the RG running factor Zy A+Av . s (/i m i n ) are a l so shown in 
Tables lC.7IG8l The errors borne by these numbers have been computed as outlined 
in Appendix B of ref. [3]. They do not include the effect of the uncertainty in 
the determination of L max /ro, reported in ref. |24| . which is numerically well below 
the error already present. The most important overall feature of these results is 
that all possible fits provide numbers compatible within la. This shows that the 
fit systematics are well under control. We conservatively take our best result to 
be that of fit D, which has the largest error, and report our final best estimates 
for Zy A+Av . s (fj> m i n ) in Tabled The result of fit D for the SSF (errors included) is 
represented in the form of shaded areas in Figs. IU"3IC.4I (left columns). 

We now return to our earlier discussion concerning the systematic effect on the 
two factors on the r.h.s. of eq. (|3.18j) . induced by the poor NLO behaviour of the 
perturbative series in some schemes. By comparing final results obtained with S2 
fixed to the NLO value to those where S2 is a free fitting parameter, we confirm that 
the use of perturbative input for the fit introduces no significant effect to the first 
factor (i.e. the running between the scsiles //mm 

and 2 /x m i n ). A proper assessment 
of the systematics on the second factor could only be obtained by calculating it 
to NNLO, which in turn would require knowledge of the perturbative coefficient 

±(2) 

7va+av; S - As the former is not available, we can estimate the size of the effect 

±(2) 

by redoing the computation with an educated guess for 7 va+ av;s- We have used 
two ansatze: First, we postulate that 7va+av; S /7va + av ;s = 7va+av ; s/7va+av- Second, 
7va+av;s is obtained from the perturbative expression 



«3 = 7vA+AV; S ( ln2 ) + [7VA+AV7VA+AV; S + 2 ^07vA+AV;s + &l7vA+Av| ( ln 2 

[1/ ±(0) n3,l/±(0) v 2 4 2 ±(0) l/i_ 9 \3 
1.6 v t 'o(7vA+AvJ + gO 7 VA+AV J (in2J 



±(0) ±(1) . „, +(1) . , ±(0) 



,2 



(4.9) 



with estimated from Fit C. The outcome of both checks is that the running factors 
in Tabled for Ova+av remain compatible within errors for all schemes. In the case of 
Ova+av they change by more than one standard deviation only for s = 1, 3, 7. These 
are indeed the schemes with largest NLO anomalous dimensions. We then conclude 
that the systematic uncertainty induced by the NLO matching in these three cases 
is not safely covered by the quoted error, and therefore these schemes should be 
discarded. As an even more conservative approach, we suggest that all schemes for 
which |7va+av;s/7va+avI > 0.2 be discarded. This stricter requirement would leave 
us with schemes 1,3,7 for Ova+av an d schemes 2,4,5,6,8,9 for Ov A+AV , which is our 
final choice of schemes deemed fully reliable. 

The RG running of the two operators is shown in Figs. IC.3IC. 41 fright columns). 
A few comments are in place: 

(i) What is plotted is the RG running of the inverse of ^vA+AV-s' which has the 
same scale dependence as the physical matrix elements of the corresponding 
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VA+AV;s(/ i min) ^vA+AV;s(^min) 



1 


1.111(19) 


0.486(7)* 


2 


1.074(24)* 


0.451(10) 


3 


1.008(19) 


0.398(7)* 


4 


1.190(24)* 


0.541(9) 


5 


1.171(23)* 


0.522(10) 


6 


1.315(24)* 


0.549(9) 


7 


1.151(19) 


0.453(7)* 


8 


1.358(25)* 


0.618(10) 


9 


1.338(23)* 


0.598(9) 



Table 1: Final results for the RG running factors Zy A+Av . s (fi m i n ). The schemes 
which suffer from systematic uncertainties related to perturbation theory have been 
indicated with an asterisk (as argued in the text, a strict criterion adopted for 
discarding a scheme is |7va+av :s /7va+av| > 0.2). 

operator (cf. eq. (|2.17|) ). 

(ii) A glance at eq. (|3.18j) reminds us that ^va+av-s(^) * s the product of the evo- 
lution function Uy A+Av . s (2 k [i m i n , ji) and the quantity ^y A+Av . s (2 fc ^ m i n ). While 
the latter quantity is computed in PT (with a 3-loop /^-function and a 2-loop 
anomalous dimension), the former is the key outcome of our non-perturbative 
calculation. Thus, by construction, the non-perturbative points coincide at 
scale 2 k fi m - m with the perturbative curve, evaluated at the same order in PT 
as the quantity ^vA+AV;s( 2fe 

Mmin ) • 

(hi) The two perturbative curves in each plot are independent of any parameters, 
once the scale and the coupling are fixed. The degree of convergence of the two 
curves at large scales n/A reflects the reliability of the perturbative estimates 
of the operator anomalous dimension. Clearly, some schemes show a better 
perturbative behaviour than others. 

(iv) The non-perturbative points are obtained as in eq. (|3.18jl . with the factor 
Uy A+Av . B (2 k /i m i n , n) calculated from fit D of the SSF. In some cases the non- 
perturbative result follows closely the NLO perturbative one up to surprisingly 
small scales. This is explicitly seen to be a scheme dependent situation. 

(v) These plots justify our strict criterion of scheme selection, as detailed above. 
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5 Connection to hadronic observables 



The RGI operator, as defined in eq. (|2.17|) . can be connected to its bare counterpart 
via a total renormalization factor, given by 

OyA+AvOr) = ^VA + AV ; s(fl'0)OvA+Av(^;50) • (5-1) 

Once the RG running of the four-fermion operator from the reference scale //min — 
max ) has been determined via the SSF, this factor decomposes into: 

-^va+av; S (3o) = Zy A+AY . s (fi min ) Z^ A+Av . s (g ,afj, m i n ) . (5.2) 

We stress that 2^ A+AV . B is a scale-independent quantity, which furthermore depends 
on the renormalization scheme only via cutoff effects. On the other hand, it de- 
pends on the particular lattice regularization chosen, though only through the factor 
Zy A+Av . s (go, aij, m i n ) , the computation of which is much less expensive than the one 
of the running Z^ A+AV;s (^ min ). 

The non-perturbative computation of Z^ A+Av . B (gQ, a/i mm ) has been performed 
at four values of (3 for each scheme and four-fermion operator, both with Clover 
and Wilson actions. The results are given in Tables IC~9llC .121 Upon multiplying by 
the corresponding ratios in Tabled the total renormalization factors are obtained. 
These can be further fitted to polynomials of the form 

^va + av; S (3o) =af + bf((3 -Q) +C f((3- 6) 2 , (5.3) 

which can be subsequently used to obtain the total renormalization factor at any 
value of (3 within the covered range [6.0219,6.4956], which comprises the typical (5- 
values used in the computation of bare observables in physically large volumes. 8 We 
supply in Table [2] the resulting fit coefficients for both the Clover and the Wilson 
case. These parameterisations represent our data with an accuracy of at least 1% 
(this comprises the point (3 = 6.0). The contribution from the error in the RGI 
renormalization factors of Tabled has been neglected: since these factors have been 
computed in the continuum limit, they should be added in quadrature after the 
quantity renormalized with the factor in eq. H5.2|) has been extrapolated itself to the 
continuum limit. 



For f3 — 6.0 a short extrapolation is necessary. 
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Clover action Wilson action 

f bf cf af bf t 



1 


0.884 


0. 


17 


0.00 


0. 


.710 


0.24 


-0.05 


2* 


0.929 


0. 


13 


0.06 


0. 


.725 


0.26 


-0.05 


3 


0.890 


0. 


15 


0.02 


0. 


.701 


0.26 


-0.07 


4* 


0.939 


0. 


14 


0.05 


0. 


.736 


0.25 


-0.03 


5* 


0.930 


0. 


14 


0.04 


0. 


.724 


0.25 


-0.01 


6* 


0.925 


0. 


15 


0.04 


0. 


.741 


0.21 


0.02 


7 


0.886 


0. 


17 


0.00 


0. 


.710 


0.23 


-0.02 


8* 


0.934 


0. 


15 


0.04 


0. 


.745 


0.22 


0.01 


9* 


0.926 


0. 


15 


0.03 


0. 


.734 


0.21 


0.03 



Clover action Wilson action 



s 


a s 


K 


c s 


a s 


K 




1* 


0.267 


0.02 


-0.03 


0.311 


-0.06 


0.05 


2 


0.293 


-0.01 


0.03 


0.324 


-0.06 


0.06 


3* 


0.269 


0.01 


-0.01 


0.308 


-0.04 


0.01 


4 


0.301 


-0.01 


0.02 


0.329 


-0.06 


0.08 


5 


0.288 


-0.01 


0.02 


0.318 


-0.06 


0.08 


6 


0.290 


-0.01 


0.02 


0.329 


-0.08 


0.09 


7* 


0.266 


0.02 


-0.02 


0.311 


-0.05 


0.03 


8 


0.299 


-0.01 


0.01 


0.334 


-0.08 


0.10 


9 


0.288 


0.00 


0.01 


0.323 


-0.08 


0.09 



Table 2: Fits to the total renormalization factor of eq. 1)5.2(1 . The schemes which suf- 
fer from systematic uncertainties related to perturbation theory have been indicated 
with an asterisk (cf. Section EJ). 



19 



6 Conclusions 



The present work is the first non-perturbative calculation of the RG evolution func- 
tion of four-fermion operators for scales ranging from the hadronic to the perturba- 
tive regime. We limit ourselves to operators with a "left-left" Dirac structure, which 
are multiplicatively renormalizable. This is the simplest possible operators 
with other Dirac structures mix under renormalization. 

The method employed is the finite size scaling approach based on the Schrodin- 
ger Functional. The Wilson lattice regularization has been used for both gluon and 
fermion fields. Combining lattice results from Wilson and Clover fermion actions 
enhances our control of continuum limit extrapolations, when obtaining the contin- 
uum step scaling function for a large range of scales. From the step scaling function 
and the perturbative estimate of the operator anomalous dimension at NLO, we 
obtain the ratio of the RGI operator to its renormalized counterpart at a hadronic 
scale. Nine different renormalization schemes have been used for each operator. 
Some of these schemes have turned out to be unstable, but this is only due to the 
bad convergence of the perturbative result for the anomalous dimension at NLO. 

We envisage that our results will be used as follows: 

(i) In simulations using Wilson type fermions and the bare operators of eq. eq. 1)3.1)1 
the matrix elements of these bare operators at fixed (3 should be multiplied 
by the renormalization factors given in eq. (|5.3[) . including a 1 percent error 
in quadrature. After continuum extrapolation, an additional error should be 
included in quadrature, corresponding to the errors quoted in tabled 

(ii) In simulations using some variant of Ginsparg- Wilson quarks (overlap quarks, 
domain-wall quarks, etc.) the results of tabled can still be used, as these are 
obtained in the continuum limit. What needs to be re-done is the calculation 
of the renormalization factor at the low energy matching scale L = 1.436 ro 
(the equivalent of tables 13 and 14). There are two ways of achieving this: 

• Via a direct evaluation of the renormalization conditions of eq. 1)3.8 J) 
at the matching scale. Obviously, this requires the formulation of the 
Schrodinger functional for Ginsparg- Wilson type quarks, e.g. along the 
lines of ref. [3*0] . 

• Via an indirect matching, as done in [3^ f° r the chiral condensate. In 
order to achieve this one just needs to compute in both regularizations a 
matrix element of the four-quark operator at matched physical conditions. 
The ratio between the bare matrix element computed with Ginsparg- 
Wilson fermions and the renormalized matrix element with Wilson quarks 
(in a given SF scheme) then yields the desired matching factor. 

A first application using Wilson type fermions consists in the computation of 
Bk in a tmQCD framework. Preliminary results have appeared in ref. 132) . 
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Appendix A A note on the difference between SF schemes 

The renormalization constants computed perturbatively at one-loop in |14j are equal 
for some of the SF schemes considered in the present work, namely for schemes s = 
1, 7 for both Ova+av an d Ova+av Consequently, the same is true for the respective 
NLO anomalous dimensions. This may suggest that the two schemes are identical. 

The non-perturbative results of this work show no such identity at the level of 
renormalization constants: at the largest values of the renormalized coupling, the 
values of Z VA+AV1)7 typically differ by several standard deviations. The difference 
is however less marked for the step scaling functions Sy A+Avl 7 , which even at the 
strongest couplings differ only by around la. 

This leaves us with the possibility that the two anomalous dimensions (which 
are defined in the continuum limit) are identical. Our data do not allow to discard 
this possibility, since the two SSFs exhibit good compatibility (cf. Tables RT5llC6|) . 
It has to be noted, however, that the continuum limit result is remarkably similar for 
many of the schemes considered. Therefore, the question whether identities between 
different schemes take place in the continuum limit cannot be strictly decided based 
on the available data. 



Appendix B Continuum limit extrapolation 

The results of the fitting procedures adopted can be summarised as follows: 

(i) The statistical accuracy of our result for 0va+av;s is always better than 2% 
and typically of 0(1 — 2%) for the largest couplings. The results for the linear 
or quadratic coefficients p have large statistical uncertainties (up to 100%), 
reflecting an overall weak cutoff dependence of Ev A+Av . s . 
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(ii) The results for <Tva+av ;s obtained by a 3-point fit are compatible to those ob- 
tained by a 4-point fit (at fixed coupling u), when the Clover action is used. 
This is also true for o\7a+av ;s > with only a few exceptions (schemes s = 3,7), 
where for one or two couplings there are discrepancies of at most 1.3a. 

With the Wilson action the situation tends to worsen. For <7va+av ;s we have 
(for each scheme) up to two or three couplings which show discrepancies of at 
most 1.6a, while for o\7a+av ; s we have up to six couplings (depending on the 
scheme) which show discrepancies typically ranging from la to 3.3a. 

We do not see any systematic trend related to the fitting ansatz (linear or 
quadratic) . 

Naturally, 3-point fit results have a larger error. 

(iii) The results for <7y A+AV;s obtained by fitting 3-points linearly are always compat- 
ible to those obtained by a quadratic fit (at fixed coupling u) when the Clover 
action is used. When 4 points are fitted, there is occasional disagreement (at 
worst for two couplings and 3a for most schemes). With the Wilson action 
things are less stable: with 4-point fits there are discrepancies for up to seven 
couplings per scheme (worst case is 6a at strong coupling). With 3-point fits 
we have up to three discrepancies per scheme (worst case is 4a). 

The results for o\7a+av;s with the Clover action show marked disagreement for 
up to nine couplings per scheme between linear and quadratic fitting (worst 
case is 6a), when 4-point fits are used. With 3-point fits we have at worst 
discrepancies at three couplings (schemes 3,7) at the 2a level. The Wilson 
results show discrepancies (up to 6a) for most couplings irrespective of the 
number of fitted points. 

(iv) The goodness of fit is satisfactory (x 2 /d.o.f. < 3) in most cases, while in a 
limited number of couplings the value tends to rise considerably. This does 
not depend systematically on the number of fitted points and choice of fitting 
ansatz. In any case, given the small number of fitted data points, x 2 /d.o.f. is a 
goodness-of-fit criterion of relatively limited value. Instead, the total x 2 /d.o.f. 
varies mostly between 1 and 2, indicating satisfactory overall quality of the fits, 
save for a few exceptions for o\7a+av; S (Wilson case with 4-point fits) where the 
value is as high as 5. 

One-loop discretisation effects can be divided out of the lattice SSF by defining 
the quantity 

Vi± l 1T\ ^VA+AV;s( U > a / L) m ^ 

£ V A+AV;>, a/ L ) = 7- g , 1 j\ ■ I 6 - 1 ) 

1 + uk^a/L) 



22 



The coefficient kf. s (a/L) is denned by expanding the ratio Sva+avjs/cva+avis i n per- 
turbation theory as 

^ A r v;s( T( L) = 1 + u k U a / L ) + - 2 k U a / L ) + ■ ■ ■ ( b - 2 ) 

and has been computed at various values of a/L in [Tl], where it is given in terms 
of the quantity Sf = kf. s /( , jy}^ Av ln2). The continuum limit of Xva+av ;s is trivially 
the same as that of £va+av ; s> but the former quantity may approach it faster, as it 
has discretisation errors which are of order u 2 . 

We find that the above procedure has significant impact on the Wilson case. The 
fits become more stable in several ways which are discussed here in correspondence 
to the criteria listed above: 

(ii) The results for CTva+av ;s obtained by a 3-point fit are incompatible to those 
obtained by a 4-point fit (at fixed coupling it) only in three schemes for at 
most 5 couplings and with a 1.5a discrepancy. 

(iii) The results for <Tva + av ;s obtained by fitting 3-points linearly show discrepancies 
to those obtained by a quadratic fit (at fixed coupling u) for at most 3 couplings 
per scheme. These discrepancies are typically 2a (in one case 3a). For o"va+av ;s 
and for two schemes only, we have discrepancies of less than 2a for a few 
couplings. 

(iv) The total x 2 /d.o.f. is always below 1.5. 

For the Clover case the continuum extrapolation of Sva+av is always compatible 
to that of Sva+av Furthermore, no significant change in the error size of the ex- 
trapolated results has been observed. For the Wilson case, where perturbative cutoff 
effects are in general large, the slope of the extrapolation decreases quite significantly 
in most cases and certainly at strong couplings. However, the extrapolated values 
from Sva+av and £va+av are again compatible and bear similar errors, but for a few 
exceptions (strongest couplings in schemes 1,3,7), where the difference between the 
extrapolated values from Xva+av an d £va+av is slightly larger than la. We conclude 
on the grounds of the above considerations, that the best result for 0va +A v;s 

(u) in 

the Wilson case is that obtained by extrapolating Sva+av 

Finally, in the spirit of ref. [5], we perform combined fits of Clover and Wilson 
data (at fixed renormalized coupling), constrained to a common continuum limit. 
This is expected to reduce the uncertainty of the results for a VA+AV . To muster 
support for this procedure we have checked the compatibility of the values of c^a_|_ av , 
obtained from linear three-point fits to the Clover data, to those obtained from linear 
three-point fits to perturbatively 0(a) improved Wilson data; recall that these are 
our best fits for each of the two datasets. In each renormalization scheme the two 
results only disagreed (typically by 1 to 2 a and at worst by 2.5 a) in a few cases 
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(for one, two or three couplings). These rare discrepancies appear both at weak and 
strong couplings. Overall, this is supportive of the universality of the continuum 
limit and justifies the option of constrained fits. For these fits the typical x 2 /d.o.f . 
range is between 1 and 2 and at worst 4, while the total x 2 /d.o.i. is around 1.2 for 
^va+av an d 1-0 for S VA+AV . 
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Table C.l: Results for the step scaling function Sv A+ av ; i with Clover action. 
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Table C.2: Results for the step scaling function Ev A +av ; i with Wilson action. 
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o 


7827 


12) 


745719) 


9527C28) 


9.7500 


8 


1.0989 


13) 


o 


131312(3) 


o 


7736 


11) 


741917) 


9590f26) 


10.0577 


12 


1.0989 


10) 


o 


131079(331 


o 


7654 


17) 


7359(21) 


9615(35) 

V7«(7WJ-t7y<JtJ / 


10.3419 


16 


1.0989 


14) 





130876(22) 





7645 


18) 


0.7320(21) 


0.9575(36) 


8.8997 


6 


i 9/1 


'i ^ 

lo J 


o 


132072(9) 


o 


7598 


8) 


7256(15) 


9550f22) 


9.1544 


8 


1.2430 


14) 


o 


131838(4) 


o 


7562 


11) 


7188(21) 

V7 • 1 1U<J \ iil / 


9505 ( 31) 


9 5202 


12 


1.2430 


;35) 


n 


131503(3) 


fl 


7473 


1 31 


7159f23) 


9580 C 35) 


9.7350 


16 


1.2430 


;34) 





131335(3 





7444 


22) 


0.7145(24) 


0.9598(43) 


8.6129 


6 




'1 9 s ! 


o 


132380(6) 


o 


7473 


13) 


7092f21) 


9490f33) 


8 8500 


8 


1.3293 


;2i) 


o 


132 140 ( 5) 


o 


7416 


13) 


7096f21) 


9569 f33) 


9.1859 


12 


1.3293 


;60) 


o 


131814(3) 


o 


7367 


19) 


6985f24) 


9481C41) 


9.4381 


16 


1.3293 


10) 





131589(2) 





7303 


17) 


0.6983(34) 


0.9562(52) 


8.3124 


6 


1.4300 


;20) 





132734(10) 





7371 


10) 


0.6939(16) 


0.9414(25) 


8.5598 


8 


1.4300 


[21) 





132453(5) 





7308 


12) 


0.6901(21) 


0.9443(33) 


8.9003 


12 


1.4300 


[50) 





132095(3) 





7219 


19) 


0.6869(27) 


0.9515(45) 


9.1415 


16 


1.4300 


[58) 





131855(3) 





7196 


15) 


0.6850(27) 


0.9519(42) 


7.9993 


6 


1.5553 


;!5) 





133118(7) 





7187 


9) 


0.6723(24) 


0.9354(35) 


8.2500 


8 


1.5553 


;24) 





132821(5) 





7135 


14) 


0.6747(24) 


0.9456(38) 


8.5985 


12 


1.5533 


70) 





132427(3) 





7109 


20) 


0.6681(29) 


0.9398(49) 


8.8323 


16 


1.5533 


70) 





132169(3) 





7023 


23) 


0.6667(34) 


0.9493(58) 



Table C.3: Results for the step scaling function S VA+AV;8 with Clover action. 
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*8~ («>>§) Zs(9o,f) Eg(u,£) 



7.7170 


6 


7.9741 


8 


8.3218 


12 


8.5479 


16 


7.4082 


6 


7.6547 


8 


7.9993 


12 


8.2415 


16 


7.1214 


6 


7.3632 


8 


7.6985 


12 


7.9560 


16 


6.7807 


6 


7.0197 


8 


7.3551 


12 


7.6101 


16 


6 551 2 




6.7860 


8 


7.1190 


12 


7.3686 


16 


6.3665 


6 


6.6100 


8 


6.9322 


12 


7.1911 


16 


6.2204 


6 


6.4527 


8 


6.7750 


12 


7.0203 


16 



1.6950(26 
1.6950(28 
1.6950(79 
1.6950(90 

1.8811(22 
1.8811(28; 
1.8811(38; 
1.8811(99 

2.1000(39 
2.1000(45; 
2.1000(80 
2.100(11) 

2.4484(37 
2.4484(45; 
2.4484(80; 
2.448(17) 

2.770(7) 
2.770(7) 
2.770(11) 
2.770(14) 

3.111(4) 
3.111(6) 
3.111(12) 
3.111(16) 

3.480(8) 
3.480(14) 
3.480(39) 
3.480(21) 



0.133517(8) 
0.133179(5) 
0.132756(4) 
0.132485(3) 

0.133961(8) 
0.133632(6) 
0.133159(4) 
0.132847(3) 

0.134423(9) 
0.134088(6) 
0.133599(4) 
0.133229(3) 

0.134994(11) 
0.134639(7) 
0.134141(5) 
0.133729(4) 

0.135327(12) 
0.135056(8) 
0.134513(5) 
0.134114(3) 

0.135488(6) 
0.135339(3) 
0.134855(3) 
0.134411(3) 

0.135470(15) 
0.135543(9) 
0.135121(5) 
0.134707(4) 



0.7045(13) 
0.6998(11) 
0.6964(7) 
0.6887(23) 

0.6828(14) 
0.6801(13) 
0.6775(8) 
0.6762(32) 

0.6622(14) 
0.6583(13) 
0.6568(8) 
0.6500(21) 

0.6330(16) 
0.6295(14) 
0.6287(9) 
0.6280(23) 

0.6079(17) 
0.6056(15) 
0.6057(10) 
0.6059(27) 

0.5830(18) 
0.5868(17) 
0.5860(22) 
0.5931(21) 

0.5615(9) 
0.5629(11) 
0.5703(20) 
0.5671(20) 



0.6558(8) 
0.6551(19) 
0.6533(25) 
0.6472(33) 

0.6307(9) 
0.6336(23) 
0.6273(30) 
0.6268(27) 

0.6036(10) 
0.6053(24) 
0.6060(27) 
0.6027(24) 

0.5639(10) 
0.5667(27) 
0.5669(29) 
0.5737(26) 

0.5311(8) 
0.5403(30) 
0.5398(21) 
0.5396(30) 

0.4976(26) 
0.5144(32) 
0.5090(29) 
0.5226(34) 

0.4652(24) 
0.4739(35) 
0.4877(31) 
0.4811(32) 



0.9309(21) 
0.9361(31) 
0.9381(37) 
0.9397(57) 

0.9237(23) 
0.9316(38) 
0.9259(46) 
0.9269(59) 

0.9115(24) 
0.9195(41) 
0.9227(43) 
0.9272(48) 

0.8908(28) 
0.9002(47) 
0.9017(48) 
0.9135(53) 

0.8737(28) 
0.8922(54) 
0.8912(38) 
0.8906(63) 

0.8535(52) 
0.8766(60) 
0.8686(59) 
0.8811(65) 

0.8285(45) 
0.8419(64) 
0.8552(62) 
0.8484(64) 



Table C.3: (continued) 
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a 


9\L) 






Z 8 (90 


'I) 


Z 8 (90, f) 


S 8 («, f) 


10 7503 


6 


U.OO l o 


'5) 


o 


134696(71 


0.8410 


11) 


789114) 


9383(21) 


11 0000 


8 


0.8873 


10) 


o 


134548(6) 


0.8248 


10) 


780915) 


9468(22) 


11.3384 


12 


0.8873 


;30) 


o 


134277(5) 


0.8110 


13) 


7709(20) 


9506(29) 


11.5736 


16 


0.8873 


;25) 





134068(6) 


0.8024 


13) 


0.7691(24) 


0.9585(34) 


10 0500 


6 


fl QQzLzL 


'7) 
J) 


o 


135659(8) 


0.8243 


12) 


7698 ( 15) 


9339(23) 


10 3000 


8 


0.9944 


13) 


o 


135457(5) 


8093 


10) 


758617) 


9374(24) 


10 6086 


12 


0.9944 


;30) 


o 


1351601) 


0.7956 


14) 


7531(26) 


9466(37) 


10.8910 


16 


0.9944 


;28) 





134849(6) 


0.7876 


15) 


0.7510(30) 


0.9535(42) 


9 5030 


6 


1 DQ8Q 




o 


136520(5) 


0.8121 


12) 


752117) 


9261(25) 


9.7500 


8 


1.0989 


13) 


o 


136310(3) 


0.7966 


11) 


743916) 


9338(24) 


10.0577 


12 


1.0989 


10) 


o 


135949(4) 


0.7818 


15) 


7341(25) 


9390(37) 


10.3419 


16 


1.0989 


14) 





135572(4) 


0.7739 


22) 


0.7293(23) 


0.9424(40) 


8.8997 


6 




'i ^ 


o 


137706(5) 


0.7933 


13) 


7321(17) 

\J • 1 '.7 _ 1111 / 


9229(26) 


9.1544 


8 


1.2430 


14) 


o 


137400(4) 


0.7786 


11) 


7188(20) 


9232(29) 


9 5202 


12 


1.2430 


;35) 


o 


136855(2) 


0.7609 


15) 


7120(22) 


9357(34) 


9.7350 


16 


1.2430 


;34) 





136523(4 


0.7514 


17) 


0.7084(28) 


0.9428(43) 


8.6129 


6 




} 


o 


138346(6) 


0.7872 


15) 


7158fl5) 


9093(26) 


8 8500 

L7 • L7 ^y V7 


8 


1.3293 


[21) 


o 


138057(4) 


0.7694 


12) 


7038(20) 


9147(30) 


9.1859 


12 


1.3293 


\m) 


o 


137503 ( 2) 


0.7503 


16) 


6980(23) 


9303(37) 


9.4381 


16 


1.3293 


10) 





137061(4) 


0.7430 


25) 


0.6891(29) 


0.9275(50) 


8.3124 


6 


1.4300 


;20) 





139128(11) 


0.7693 


15) 


0.7017(18) 


0.9121(29) 


8.5598 


8 


1.4300 


[21) 





138742(7) 


0.7569 


13) 


0.6885(21) 


0.9096(32) 


8.9003 


12 


1.4300 


;50) 





138120(8) 


0.7419 


18) 


0.6877(32) 


0.9269(49) 


9.1415 


16 


1.4300 


;58) 





137655(5) 


0.7316 


19) 


0.6818(31) 


0.9319(49) 


7.9993 


6 


1.5553 


;!5) 





140003(11) 


0.7600 


16) 


0.6797(19) 


0.8943(31) 


8.2500 


8 


1.5553 


;24) 





139588(8) 


0.7423 


13) 


0.6739(22) 


0.9079(34) 


8.5985 


12 


1.5533 


70) 





138847(6) 


0.7300 


19) 


0.6681(35) 


0.9152(54) 


8.8323 


16 


1.5533 


70) 





138339(7) 


0.7208 


28) 


0.6594(31) 


0.9148(56) 



Table C.4: Results for the step scaling function E VA+AV;8 with Wilson action. 
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P 


a 


9 (L) 




r, / T,\ 


Z 8 {90, f) 


S 8 («, f ) 


7.7170 


6 


1 6950(261 


140954(12) 


7443(16) 


6638(20) 


8918(33) 


7.9741 


8 


1 6950(28) 


140438(8) 


7292(15) 


6636(30) 


9100(45) 

VX • »_/ 1 V7 VX 1 a ^X / 


8.3218 


12 


1 6950(79) 


139589(6) 


7131(20) 


6535(30) 


9164(49) 


8.5479 


16 


1.6950(90) 


0.139058(6) 


0.7038(22) 


0.6445(41) 


0.9157(65) 


7.4082 


6 


1 8811(22) 


142145(11) 


7280(17) 


6404(21) 


8797(35) 


7.6547 


8 


1 8811(28) 

i- . W \-J 1 1 \ vX 1 


141572(9) 


7123(16) 


6339(25) 


8899(40) 


7.9993 


12 


1 8811(38) 


140597(6) 


6965(21) 


6269(33) 


9001(55) 


8.2415 


16 


1.8811(99) 


0.139900(6) 


0.6855(26) 


0.6303(36) 


0.9195(63) 


7.1214 


6 


2 1000(39) 


143416(11) 


7040(18) 


6116(23) 


8688(40) 

vy • w Vy W LJ 1 A \J 1 


7.3632 


8 


2 1000(45) 


142749(9) 


6886(16) 


6009(32) 


8726(51) 

KJ *\J 1 ^Ul *J1 / 


7.6985 


12 


2 1000(80) 


141657(6) 


6754(16) 


5992(22) 


8872(39) 


7.9560 


16 


2.100(11) 


0.140817(7) 


0.6688(24) 


0.5995(29) 


0.8964(54) 


6.7807 


6 


2 4484(37) 


145286(11) 


6818(18) 

U . UU 1 U 1 1 LJ / 


5807(22) 


8517(39) 


7.0197 


8 


2 4484(45) 


n 144454(7) 

V / • -L A l i '..7 1 i | / 


6644(17) 


5678(25) 


8546(44) 


7.3551 


12 


2 4484(80) 


143113(6) 


6537(19) 

\j • xj^jtj i y it/ / 


5576(35) 


8530(59) 


7.6101 


16 


2.448(17) 


0.142107(6) 


0.6452(25) 


0.5735(42) 


0.8889(74) 


6.5512 


6 


2.770(7) 


146825(11) 


6596(21) 

U>V/Uf/U \ ill / 


5459(23) 


8276(44) 


6.7860 


8 


2.770(7) 


145859(7) 


6430(20) 


5370(34) 


8351(59) 


7.1190 


12 


2 770(11) 


144299(8) 

vy • 1 A A^ I U / 


6303(24) 


5366(28) 


8513(55) 


7.3686 


16 


2.770(14) 


0.143175(7) 


0.6273(31) 


0.5423(38) 


0.8645(74) 


6.3665 


6 


3.111(4) 


0.148317(10) 


0.6420(22) 


0.5124(21) 


0.7981(43) 


6.6100 


8 


3.111(6) 


0.147112(7) 


0.6276(21) 


0.5100(25) 


0.8126(48) 


6.9322 


12 


3.111(12) 


0.145371(7) 


0.6124(22) 


0.5053(38) 


0.8251(69) 


7.1911 


16 


3.111(16) 


0.144060(8) 


0.6059(27) 


0.5151(33) 


0.8501(66) 


6.2204 


6 


3.480(8) 


0.149685(15) 


0.6236(24) 


0.4871(25) 


0.7811(50) 


6.4527 


8 


3.480(14) 


0.148391(9) 


0.6017(22) 


0.4759(43) 


0.7909(77) 


6.7750 


12 


3.480(39) 


0.146408(7) 


0.5955(21) 


0.4808(46) 


0.8074(82) 


7.0203 


16 


3.480(21) 


0.145025(8) 


0.5863(25) 


0.4917(45) 


0.8386(85) 



Table C.4: (continued) 
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Figure C.l: Continuum extrapolations of Ev A+A v ; i a t fixed renormalized coupling u 
for the improved action (full symbols, solid line) and the unimproved action (open 
symbols, dashed line). The L/a = 6 data points have not been included in the fits. 
The value of u increases from top to bottom and from left to right. 
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Figure C.2: Continuum extrapolations of £y A+AV;g a ^ fixed renormalized coupling u 
for the improved action (full symbols, solid line) and the unimproved action (open 
symbols, dashed line). The L/a = 6 data points have not been included in the fits. 
The value of u increases from top to bottom and from left to right. 
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u 


at(u) 


a+(u) 


a£(u) 


ai{u) 


°t(u) 


0.8873 


1.024(6) 


1.022(7) 


1.026(7) 


1.020(7) 


1.020(7) 


0.9944 


1.026(7) 


1.028(8) 


1.029(7) 


1.025(8) 


1.025(8) 


1.0989 


1.011(7) 


1.008(8) 


1.012(7) 


1.006(7) 


1.007(7) 


1.2430 


1.030(8) 


1.031(9) 


1.034(8) 


1.028(9) 


1.028(8) 


1.3293 


1.030(8) 


1.031(10) 


1.034(9) 


1.027(9) 


1.027(9) 


1.4300 


1.031(9) 


1.035(11) 


1.035(9) 


1.031(10) 


1.031(10) 


1.5553 


1.034(10) 


1.038(12) 


1.040(10) 


1.032(11) 


1.033(11) 


i.oyou 


i.u4o^y j 




1.1)00^111 J 


l.U4o(ll) J 


1 C\A A ( ~\ C\\ 

l.U44^1U J 


1.8811 


1.030(10) 


1.035(12) 


1.039(10) 


1.027(11) 


1.028(11) 


2.1000 


1.045(9) 


1.060(12) 


1.058(10) 


1.048(10) 


1.050(10) 


2.4484 


1.074(10) 


1.093(13) 


1.095(11) 


1.073(12) 


1.074(12) 


2.770 


1.104(11) 


1.143(15) 


1.138(12) 


1.110(13) 


1.112(12) 


3.111 


1.118(12) 


1.157(17) 


1.157(14) 


1.120(14) 


1.124(14) 


3.480 


1.130(12) 


1.186(18) 


1.180(15) 


1.140(15) 


1.146(14) 



u 


a£(u) 


a+(u) 




a+(u) 


0.8873 


1.018(7) 


1.023(6) 


1.017(7) 


1.017(7) 


0.9944 


1.022(8) 


1.025(7) 


1.021(8) 


1.021(8) 


1.0989 


1.005(7) 


1.010(7) 


1.004(7) 


1.004(7) 


1.2430 


1.024(8) 


1.029(8) 


1.023(8) 


1.024(8) 


1.3293 


1.024(9) 


1.029(8) 


1.023(9) 


1.022(8) 


1.4300 


1.025(9) 


1.029(8) 


1.025(9) 


1.024(9) 


1.5553 


1.027(10) 


1.033(9) 


1.025(10) 


1.026(10) 


1.6950 


1.035(10) 


1.044(9) 


1.031(10) 


1.033(9) 


1.8811 


1.018(10) 


1.028(9) 


1.016(10) 


1.017(10) 


2.1000 


1.034(9) 


1.041(9) 


1.031(9) 


1.033(9) 


2.4484 


1.050(11) 


1.066(9) 


1.045(10) 


1.047(10) 


2.770 


1.072(11) 


1.090(10) 


1.063(11) 


1.065(10) 


3.111 


1.073(12) 


1.100(11) 


1.064(11) 


1.068(11) 


3.480 


1.084(12) 


1.111(11) 


1.074(12) 


1.079(11) 



Table C.5: Continuum extrapolations of Sva+av ;s combining Clover and (pertur- 
batively O(a) improved) Wilson data. Linear dependence on (a/L) is assumed for 
both actions. The L/a = 6 data have not been taken into account. 
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u 




a 2 (u) 


cr 3 (u) 


cr 4 (u) 


°5 ( u ) 


0.8873 


0.981(5) 


0.978(6) 


0.985(6) 


0.973(5) 


0.974(6) 


0.9944 


0.971(6) 


0.975(7) 


0.976(7) 


0.971(6) 


0.970(6) 


1.0989 


0.968(5) 


0.963(7) 


0.972(7) 


0.959(6) 


0.961(6) 


1.2430 


0.976(6) 


0.980(8) 


0.984(7) 


0.974(6) 


0.975(7) 


1.3293 


0.961(6) 


0.960(9) 


0.967(8) 


0.955(7) 


0.955(8) 


1.4300 


0.963(6) 


0.975(8) 


0.972(8) 


0.968(7) 


0.969(8) 


1.5553 


0.950(8) 


0.952(11) 


0.960(10) 


0.944(8) 


0.944(9) 


1.6950 


0.954(7) 


0.962(10) 


0.970(9) 


0.949(8) 


0.951(9) 


1.8811 


0.953(7) 


0.956(11) 


0.966(9) 


0.943(9) 


0.942(10) 


2.1000 


0.934(7) 


0.967(11) 


0.959(8) 


0.947(8) 


0.951(9) 


2.4484 


0.943(8) 


0.971(13) 


0.982(10) 


0.941(10) 


0.942(10) 


2.770 


0.911(8) 


0.986(14) 


0.974(10) 


0.934(10) 


0.938(11) 


3.111 


0.908(8) 


0.987(15) 


0.982(11) 


0.930(11) 


0.939(12) 


3.480 


0.896(10) 


0.995(16) 


0.983(12) 


0.926(12) 


0.935(13) 



u 


a 6 (u) 


a 7 (u) 


(7g (u) 


a g (u) 


0.8873 


0.974(5) 


0.982(6) 


0.971(5) 


0.971(5) 


0.9944 


0.971(6) 


0.973(6) 


0.968(5) 


0.967(6) 


1.0989 


0.960(5) 


0.970(6) 


0.957(5) 


0.958(6) 


1.2430 


0.973(6) 


0.980(6) 


0.969(6) 


0.970(6) 


1.3293 


0.954(7) 


0.963(7) 


0.951(6) 


0.950(7) 


1.4300 


0.967(7) 


0.967(7) 


0.962(6) 


0.963(7) 


1.5553 


0.943(9) 


0.953(8) 


0.938(7) 


0.938(8) 


1.6950 


0.947(8) 


0.960(7) 


0.939(7) 


0.940(8) 


1.8811 


0.940(9) 


0.956(8) 


0.932(8) 


0.932(9) 


2.1000 


0.944(8) 


0.944(7) 


0.931(7) 


0.935(8) 


2.4484 


0.933(10) 


0.956(8) 


0.916(8) 


0.917(9) 


2.770 


0.925(10) 


0.932(8) 


0.895(9) 


0.898(10) 


3.111 


0.917(11) 


0.935(9) 


0.885(9) 


0.893(10) 


3.480 


0.910(12) 


0.927(10) 


0.872(11) 


0.880(12) 



Table C.6: Continuum extrapolations of Sva+av :s combining Clover and (pertur- 
batively 0(a) improved) Wilson data. Linear dependence on (a/L) is assumed for 
both actions. The L/a = 6 data have not been taken into account. 
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s 


Fit 


si 


S2 


S3 


s 4 


x 2 

d.o.f. 


^it (/^min) 


1 


A 


0.01755762 


0.00529(56) 






1.27 


1.095(16) 




B 


0.01755762 


0.0006(23) 


0.00167(81) 




1.02 


1.108(18) 




C 


0.01755762 


0.00137(2) 


0.00142(19) 




0.95 


1.105(15) 




Pi 


U.Ul / OD loz 




n nn/i 1 (kq\ 
u.uu4i^oy j 


— U.UUUO^IZ J 


1.1U 


1 111/101 




E 


0.01755762 


0.00137(2) 


0.0012(12) 


0.00006(38) 


1.03 


1.106(16) 


2 


A 


0.01755762 


0.00873(76) 






1.75 


1.050(21) 




B 


0.01755762 


-0.0006(29) 


0.0035(11) 




0.99 


1.067(22) 




C 


0.01755762 


-0.00320(2) 


0.00438(27) 




0.97 


1.076(20) 




L) 


U.Ul / 00 ( DZ 


— U.UUoo^o f ) 


u.uuyi^ to) 


— U.UUll^lO J 


1 no 
l.UZ 


l.U / 4(^z4 ) 




E 


0.01755762 


-0.00320(2) 


0.0060(16) 


-0.00055(51) 


0.95 


1.069(21) 


3 


A 


0.01755762 


0.00884(63) 






1.83 


0.990(17) 




B 


0.01755762 


0.0014(25) 


0.00275(90) 




1.22 


1.004(18) 




C 


0.01755762 


0.00218(2) 


0.00248(23) 




1.13 


1.001(15) 




Pi 


U.Ul / 00 ( DZ 


— v.\JUov{ (v) 


U.UUdq^qo ) 


— U.UUUo^lo J 


1 on 
i.zy 


i.uuo^iy ) 




E 


0.01755762 


0.00218(2) 


0.0024(13) 


0.00004(43) 


1.22 


1.001(16) 


4 


A 


0.01755762 


0.00537(66) 






1.57 


1.163(20) 




B 


0.01755762 


-0.0023(26) 


0.00281(94) 




0.95 


1.182(22) 




C 


0.01755762 


-0.00451(2) 


0.00356(23) 




0.93 


1.191(19) 




D 


01 755762 


—0 0086f80"> 


0084ffi8l 


—0 001 1 f14l 


97 


1 190(24"! 




E 


0.01755762 


-0.00451(2) 


0.0050(14) 


-0.00047(45) 


0.91 


1.183(20) 


5 


A 


0.01755762 


0.00584(63) 






1.60 


1.143(19) 




B 


0.01755762 


-0.0020(26) 


0.00285(91) 




0.92 


1.164(21) 




C 


0.01755762 


-0.00404(1) 


0.00353(22) 




0.89 


1.173(18) 




D 


0.01755762 


-0.0084(78) 


0.0085(66) 


-0.0011(13) 


0.93 


1.171(23) 




E 


0.01755762 


-0.00404(1) 


0.0049(13) 


-0.00044(43) 


0.88 


1.165(19) 



Table C.7: Fits to the continuum step scaling functions cr^ A+Av . s and results for the 
ratio Zy A+AV;s (/x m i n ) . 
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2 A. 

S Fit Sl S 2 S 3 S 4 gj^-j- Z s + (/x min ) 



6 


A 


0.01755762 


0.00136(56) 






0.96 


1.294(20) 




B 


0.01755762 


-0.0030(24) 


0.00154(83) 




0.75 


1.309(22) 




C 


0.01755762 


-0.00442(2) 


0.00203(19) 




0.72 


1.317(18) 




D 


0.01755762 


-0.0071(74) 


0.0052(62) 


-0.0007(12) 


0.79 


1.315(24) 




E 


0.01755762 


r\ r\r\ A 

—0.00442(2) 


0.0030(12) 


—0.00030(39) 


0.74 


1.310(20) 


7 


A 


0.01755762 


0.00387(51) 






1.06 


1.138(16) 




B 


0.01755762 


0.0004(22) 


0.00122(77) 




0.93 


1.148(17) 




C 


0.01755762 


0.00137(2) 


0.00091(18) 




0.87 


1.144(14) 




D 


0.01755762 


-0.0019(68) 


0.0033(56) 


-0.0004(11) 


1.00 


1.151(19) 




E 


0.01755762 


0.00137(2) 


0.0006(11) 


0.00010(36) 


0.94 


1.146(16) 


8 


A 


0.01755762 


0.00049(55) 






0.87 


1.338(20) 




B 


0.01755762 


-0.0037(24) 


0.00147(83) 




0.67 


1.353(22) 




C 


0.01755762 


-0.00533(1) 


0.00204(19) 




0.66 


1.362(18) 




D 


0.01755762 


-0.0070(73) 


0.0044(61) 


-0.0006(12) 


0.71 


1.358(25) 




E 


0.01755762 


-0.00533(1) 


0.0030(12) 


-0.00032(39) 


0.66 


1.355(20) 


9 


A 


0.01755762 


0.00092(52) 






0.92 


1.316(19) 




B 


0.01755762 


-0.0033(23) 


0.00150(79) 




0.70 


1.332(21) 




C 


0.01755762 


-0.00485(2) 


0.00200(18) 




0.67 


1.341(17) 




D 


0.01755762 


-0.0075(71) 


0.0050(59) 


-0.0007(11) 


0.72 


1.338(23) 




E 


0.01755762 


-0.00485(2) 


0.0029(12) 


-0.00030(37) 


0.68 


1.333(19) 



Table C.7: (continued) 
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s Fit 


si 


S2 


«3 


s 4 


x 2 

d.o.f. 


Z s (^min) 


1 A 


-0.03511524 


0.00252(42) 






1.87 


0.497(6) 


B 


-0.03511524 


0.0085(18) 


-0.00217(63) 




1.03 


0.490(7) 


C 


-0.03511524 


0.01531(2) 


-0.00448(15) 




2.06 


0.479(6) 


D 


-0.03511524 


0.0154(54) 


-0.0082(46) 


0.00122(91) 


0.96 


0.486(7) 


E 


-0.03511524 


0.01531(2) 


-0.00815(92) 


0.00120(30) 


0.88 


0.486(6) 



2 A 


-0.03511524 


0.00941(69) 






0.72 


0.449(9) 


B 


-0.03511524 


0.0069(26) 


0.00094(93) 




0.70 


0.450(9) 


C 


-0.03511524 


0.00539(1) 


0.00146(25) 




0.67 


0.453(8) 


D 


-0.03511524 


0.0049(77) 


0.0028(67) 


-0.0004(14) 


0.76 


0.451(10) 


E 


-0.03511524 


0.00539(1) 


0.0023(14) 


-0.00029(46) 


0.69 


0.451(9) 



3 A 


-0.03511524 


0.00934(53) 






0.75 


0.401(6) 


B 


-0.03511524 


0.0113(22) 


-0.00071(78) 




0.75 


0.399(6) 


C 


-0.03511524 


0.01612(2) 


-0.00235(19) 




1.05 


0.392(5) 


D 


-0.03511524 


0.0140(68) 


-0.0030(57) 


0.0005(11) 


0.80 


0.398(7) 


E 


-0.03511524 


0.01612(2) 


-0.0048(11) 


0.00080(37) 


0.74 


0.397(6) 



4 A 


-0.03511524 


0.00407(52) 






0.81 


0.541(8) 


B 


-0.03511524 


0.0042(20) 


-0.00006(73) 




0.87 


0.541(9) 


C 


-0.03511524 


0.00276(2) 


0.00045(18) 




0.85 


0.544(8) 


D 


-0.03511524 


0.0046(61) 


-0.0004(52) 


0.0001(11) 


0.95 


0.541(9) 


E 


-0.03511524 


0.00276(2) 


0.0012(11) 


-0.00023(35) 


0.88 


0.542(8) 



5 A 


-0.03511524 


0.00466(56) 






0.69 


0.521(9) 


B 


-0.03511524 


0.0041(23) 


0.00022(80) 




0.74 


0.522(9) 


C 


-0.03511524 


0.00451(1) 


0.00007(20) 




0.69 


0.521(8) 


D 


-0.03511524 


0.0044(67) 


-0.0001(57) 


0.0001(11) 


0.81 


0.522(10) 


E 


-0.03511524 


0.00451(1) 


-0.0002(12) 


0.00008(39) 


0.74 


0.522(8) 



Table C.8: Fits to the continuum step scaling functions cr VA+AV;s and results for the 
ratio ^ V A+AV;s(/ Xmin )- 
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s 


Fit 


si 


S2 


S3 


s 4 


x 2 

d.o.f. 


(Mmin) 


6 


A 
B 
C 
D 
E 


-0.03511524 
-0.03511524 
-0.03511524 
-0.03511524 
-0.03511524 


0.00297(52) 

0.0044(20) 

0.00416(1) 

0.0058(60) 

0.00416(1) 


-0.00052(72) 
-0.00044(18) 
-0.0018(52) 
-0.0004(11) 


0.0003(10) 
-0.00001(35) 


0.81 
0.83 
0.77 
0.90 
0.83 


0.551(9) 
0.550(9) 
0.550(8) 
0.549(9) 
0.550(8) 


7 


A 
B 
C 
D 
E 


-0.03511524 
-0.03511524 
-0.03511524 
-0.03511524 
-0.03511524 


0.00491(44) 

0.0094(19) 

0.01531(2) 

0.0154(58) 

0.01531(2) 


-0.00162(66) 
-0.00362(15) 
-0.0068(48) 
-0.00673(95) 


0.00103(95) 
0.00102(31) 


1.34 
0.95 
1.62 
0.92 
0.85 


0.462(6) 
0.457(6) 
0.447(5) 
0.453(7) 
0.453(6) 


8 


A 
B 
C 
D 
E 


-0.03511524 
-0.03511524 
-0.03511524 
-0.03511524 
-0.03511524 


-0.00017(45) 
0.0021(18) 
0.00194(1) 
0.0048(54) 
0.00194(1) 


-0.00082(64) 
-0.00078(16) 
-0.0033(46) 
-0.00089(95) 


0.00051(94) 
0.00004(31) 


1.00 
0.95 
0.88 
1.01 
0.95 


0.622(9) 

0.620(9) 

0.621(8) 

0.618(10) 

0.621(8) 


9 


A 
B 
C 
D 
E 


-0.03511524 
-0.03511524 
-0.03511524 
-0.03511524 
-0.03511524 


0.00038(50) 

0.0022(20) 

0.00370(2) 

0.0057(61) 

0.00370(2) 


-0.00065(72) 
-0.00118(18) 
-0.0038(52) 
-0.0022(11) 


0.0006(10) 
0.00032(35) 


0.89 
0.89 
0.87 
0.94 
0.87 


0.600(9) 

0.599(10) 

0.595(9) 

0.596(10) 

0.598(9) 



Table C.8: (continued) 
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Figure C.3: Left column: The step scaling function o"va+av ;s ( , u) (discrete points) as 
obtained non-perturbatively from combined fits to Clover and Wilson data. The 
shaded area is the result of fit D to the points (see text). The dotted (dashed) line is 
the LO (NLO) perturbative result. Right column: RG running of Oy A +Av obtained 
non-perturbatively (discrete points) at specific values of the renormalization scale fj,, 
in units of A (taken from ref. [I]). The lines are perturbative results at the indicated 
order for the Callan-Symanzik /3-function and the operator anomalous dimension 7. 
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Figure C.3: (continued) 
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Figure C.3: (continued) 
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Figure C.4: Left column: The step scaling function (Jy A+Av . E (u) (discrete points) as 
obtained non-perturbatively from combined fits to Clover and Wilson data. The 
shaded area is the result of fit D to the points (see text). The dotted (dashed) line is 
the LO (NLO) perturbative result. Right column: RG running of Ova+av obtained 
non-perturbatively (discrete points) at specific values of the renormalization scale /i, 
in units of A (taken from ref. [I]). The lines are perturbative results at the indicated 
order for the Callan-Symanzik /3-function and the operator anomalous dimension 7. 
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Figure C.4: (continued) 
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Figure C.4: (continued) 



46 



Clover action Wilson action 

K cr -ZvA + AV;s K cr ^VA+AV;s 
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10 


6. 


.2885 


12 


6. 


.4956 


16 



0.135043(17) 
0.135643(11) 
0.135739(13) 
0.135577(7) 

0.135043(17) 
0.135643(11) 
0.135739(13) 
0.135577(7) 

0.135043(17) 
0.135643(11) 
0.135739(13) 
0.135577(7) 

0.135043(17) 
0.135643(11) 
0.135739(13) 
0.135577(7) 

0.135043(17) 
0.135643(11) 
0.135739(13) 
0.135577(7) 



0.7985(27) 
0.8228(17) 
0.8400(20) 
0.8758(46) 

0.8661(42) 
0.8873(26) 
0.9029(31) 
0.9412(74) 

0.8851(36) 
0.9105(22) 
0.9272(26) 
0.9667(61) 

0.7911(32) 
0.8100(20) 
0.8252(24) 
0.8585(56) 

0.7963(32) 
0.8153(20) 
0.8307(24) 
0.8633(56) 



0.153371(10) 
0.152012(7) 
0.150752(10) 
0.148876(13) 

0.153371(10) 
0.152012(7) 
0.150752(10) 
0.148876(13) 

0.153371(10) 
0.152012(7) 
0.150752(10) 
0.148876(13) 

0.153371(10) 
0.152012(7) 
0.150752(10) 
0.148876(13) 

0.153371(10) 
0.152012(7) 
0.150752(10) 
0.148876(13) 



0.6435(28) 
0.6743(24) 
0.6978(33) 
0.7377(48) 

0.6802(37) 
0.7145(33) 
0.7413(46) 
0.7861(66) 

0.7010(33) 
0.7370(29) 
0.7647(40) 
0.8096(58) 

0.6232(29) 
0.6534(26) 
0.6768(36) 
0.7183(52) 

0.6226(29) 
0.6534(25) 
0.6770(35) 
0.7215(50) 



Table C.9: Results for Zy A+Av . s (go, L/a) at fixed scale L = 1.436 tq (corresponding to 
Mmin = (2^max) -1 )- Each block contains the results from a different renormalization 
condition (s = 1, . . . , 5). 
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Clover action Wilson action 

:r ^VA+AV;s K cr ^VA+AV;s 
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1628 
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.2885 
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6. 
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.2885 
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.4956 


16 


6. 


0219 
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.1628 
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6. 


.2885 


12 


6. 


.4956 


16 


6. 


0219 


8 


6. 


1628 


10 


6. 


.2885 


12 


6. 


.4956 


16 



0.135043(17) 
0.135643(11) 
0.135739(13) 
0.135577(7) 

0.135043(17) 
0.135643(11) 
0.135739(13) 
0.135577(7) 

0.135043(17) 
0.135643(11) 
0.135739(13) 
0.135577(7) 

0.135043(17) 
0.135643(11) 
0.135739(13) 
0.135577(7) 



0.7058(25) 
0.7236(16) 
0.7388(19) 
0.7686(43) 

0.7722(24) 
0.7948(15) 
0.8111(18) 
0.8446(40) 

0.6902(24) 
0.7071(15) 
0.7219(17) 
0.7501(40) 

0.6947(23) 
0.7117(14) 
0.7267(17) 
0.7543(39) 



0.153371(10) 
0.152012(7) 
0.150752(10) 
0.148876(13) 

0.153371(10) 
0.152012(7) 
0.150752(10) 
0.148876(13) 

0.153371(10) 
0.152012(7) 
0.150752(10) 
0.148876(13) 

0.153371(10) 
0.152012(7) 
0.150752(10) 
0.148876(13) 



0.5670(25) 
0.5919(22) 
0.6104(29) 
0.6486(42) 

0.6209(26) 
0.6501(22) 
0.6719(29) 
0.7122(42) 

0.5520(23) 
0.5764(20) 
0.5946(27) 
0.6319(40) 

0.5515(22) 
0.5764(19) 
0.5948(26) 
0.6347(37) 



Table C.10: Results for Zva+avjsG/Oj L/a) at fixed scale L = 1.436 ro (corresponding 
to /Lt m in = (2L max ) _1 ). Each block contains the results from a different renormaliza- 
tion condition (s = 6, . . . , 9). 
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Clover action Wilson action 



p 


L 

a 




■^VA+AV;s 




■^VA+AV;s 


6.0219 


8 


0.135043(17) 


0.5499(18) 


0.153371(10) 


0.6368(27) 


6.1628 


10 


0.135643(11) 


0.5561(11) 


0.152012(7) 


0.6212(21) 


6.2885 


12 


0.135739(13) 


0.5573(13) 


0.150752(10) 


1^1 fi -t -t V t C\C\\ 

0.6115(29) 
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U. lOOO i i \ t J 




D 1 A8£7fi/'1 ^ 

U. 14:00 / U I lOJ 


U.UUUU^U ) 
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r\ 1 OTA JO/1 ^ 

0.135043(17) 


0.6474(36) 


0.153371(10) 


0.7150(40) 


6.1628 


10 


0.135643(11) 


/"* A ^7A1 / . "\ -1 \ 

0.6472(21) 


0.152012(7) 


0.7023(32) 


6.2885 


12 


0.135739(13) 


0.6452(25) 


r\ -i r Af-7r / 1 n\ 

0.150752(10) 


0.6922(44) 


6.4956 


16 


0.135577(7) 


0.6524(57) 


0.148876(13) 


0.6895(67) 


6.0219 


8 


0.135043(17) 


0.6743(25) 


0.153371(10) 


0.7719(35) 


6.1628 


10 


0.135643(11) 


0.6814(16) 


0.152012(7) 


0.7570(28) 


6.2885 


12 


0.135739(13) 


0.6809(18) 


0.150752(10) 


0.7497(38) 


6.4956 


16 


0.135577(7) 


0.6866(41) 


0.148876(13) 


0.7311(51) 


6.0219 


8 


0.135043(17) 


0.5549(24) 


0.153371(10) 


0.6056(26) 


6.1628 


10 


0.135643(11) 


0.5538(14) 


0.152012(7) 


0.5944(21) 


6.2885 


12 


0.135739(13) 


0.5529(17) 


0.150752(10) 


0.5854(29) 


6.4956 


16 


0.135577(7) 


0.5558(37) 


0.148876(13) 


0.5869(46) 


6.0219 


8 


0.135043(17) 


0.5516(26) 


0.153371(10) 


0.6057(29) 


6.1628 


10 


0.135643(11) 


0.5514(15) 


0.152012(7) 


0.5941(23) 


6.2885 


12 


0.135739(13) 


0.5520(18) 


0.150752(10) 


0.5858(32) 


6.4956 


16 


0.135577(7) 


0.5562(40) 


0.148876(13) 


0.5868(50) 



Table C.ll: Results for Z^ A+Av . B (gQ, L/a) at fixed scale L = 1.436 vq (corresponding 
to /Lt m in = (2L max ) _1 ). Each block contains the results from a different renormaliza- 
tion condition (s = 1, . . . , 5). 
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Clover action Wilson action 



p 


L 

a 




■^VA+AV;s 




■^VA+AV;s 


6.0219 


8 


0.135043(17) 


0.5276(23) 


0.153371(10) 


0.5960(27) 


6.1628 


10 


0.135643(11) 


0.5278(14) 


0.152012(7) 


0.5819(22) 


6.2885 


12 


0.135739(13) 


0.5279(16) 


0.150752(10) 


0.5700(29) 


6.4956 


16 


0.135577(7) 


0.5328(36) 


0.148876(13) 


0.5689(44) 


6.0219 


8 


0.135043(17) 


0.5883(19) 


0.153371(10) 


0.6837(30) 


6.1628 


10 


0.135643(11) 


0.5948(12) 


0.152012(7) 


0.6677(23) 


6.2885 


12 


0.135739(13) 


0.5956(14) 


0.150752(10) 


0.6586(32) 


6.4956 


16 


0.135577(7) 


0.5999(32) 


0.148876(13) 


0.6432(41) 


6.0219 


8 


0.135043(17) 


0.4841(20) 


0.153371(10) 


0.5364(22) 


6.1628 


10 


0.135643(11) 


0.4835(12) 


0.152012(7) 


0.5243(17) 


6.2885 


12 


0.135739(13) 


0.4837(14) 


0.150752(10) 


0.5143(24) 


6.4956 


16 


0.135577(7) 


0.4856(30) 


0.148876(13) 


0.5163(38) 


6.0219 


8 


0.135043(17) 


0.4812(21) 


0.153371(10) 


0.5365(23) 


6.1628 


10 


0.135643(11) 


0.4813(12) 


0.152012(7) 


0.5241(18) 


6.2885 


12 


0.135739(13) 


0.4829(14) 


0.150752(10) 


0.5147(26) 


6.4956 


16 


0.135577(7) 


0.4860(31) 


0.148876(13) 


0.5162(40) 



Table C.12: Results for ^v A+AV;s (<7o, L/a) at fixed scale L = 1.436 ro (corresponding 
to /Lt m in = (2L max ) _1 ). Each block contains the results from a different renormaliza- 
tion condition (s = 6, . . . , 9). 
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